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The mode shapes^ associated natural frequencies 
CoUd nature and amount of damping govern the dynamic 
:/:esponse of structures, Iliis worP presents a minimiz- 
ation method of extracting the eigenvalues and eigen- 
vectors of a structural system having finite degrees of 
freedom. The structures considered ore supposed to he 
linear and undamped. The Rayleigh quotient has been 
minimized successively in reduced subspaces to obtain 
the second and higher frequencies and the mode shapes 
associated with thevii. Gonjugate Gradient Method of 
minimization has been used- Its modest demands on 
Gomputer storage space have been farther supplemented 
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by xjerfoming the analysis which does not reqaire the 
assembly of structural stiffness and mass matrices. Ihe 
method seems to be particularly useful for large sized 
problems wt'ere only first few eigenvalues and eigenvectors 
are desired, 

i'he study of free vibrations is usually an 
inteiTiiediate step to perform the forced vibration analysis, 
'Ihe method has been applied to simple structures like 
c until ever beam, portal frame and skev; bridge and 
results obtained compare very well. 



chapter - I 


IHTROPHOTIOH 


Performing a discrete element analysis of the 
structure, the matrix formulation of the general dynainic 
resi^onse problem for the undamped structure is'» 

[m] ? + [ic] Y = i’ (t) ... (1) 

where ■the n x n matrices [wi] and Ck] are respectively 
the masfj matrix and the stiffness matrix of the structure. 
These matrices are real and symmetric and their order* 
n corresponds to the degrees of freedora to which an infinite 
degree of freedom continuous structure has been discretized. 
The n X 1 vectors X , Y and P (t) represent the absolute 
acceleration, absolute displacement and load respectively. 
The set of equations Eqn. (l), consist of a set of coupled 
second order ordinary differential equations. These are 
usually solved by uncoupling the equations through matrix 
transformation. The uncoupled equations resemble the 
equation of motion of a single degree of freedom system. 

The dynomic response is finally obtained by the linear 
superp osition of the solution of these uncoupled equations. 
The matrix which performs uncoupling of the set of ‘ii)qn.(l) 
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requires ihe free vibration analysis* of the structure 
v/hich in turn means the solution to the set of equations: 

fmj Y + jic'; Y = 0 ... (2) 

^he above set of equations are reduced to an 
eigenvalue problem ' 




M I I 


. ( 3 ) 


v/‘here /\ is the sc, alar quantity called eigenvalue which 
represents the square of the natural frequency of the system 
and 1 is the (n x l) eigenvector called the principal 
mode shape associated with the particular value ^ . 


t'his work presents a method to obtain the eigen- 
values A and associated eigenvectors 'the method 

presented here is iterative and exploits the minimization 
technique o.nd is expected to be specifically useful for 
determination of the first few eigenvalues and associated 
eigenvectors for large sized structures i^e, where n is 
large* 

1,1 iiETHObS OE EIGENSOiaflOH: 

by eigensolution is meant the solution of the 
eigenproblem, Eqn. (3), which in turn means the deter- 
mination of the scalar quantities X and their associated 
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m odes hap es 

Various methods of eigonsolution known can be 
Categorized as: 

(a) '•^ransfomation I'Aethods 

( b) Iterative Methods, 

Iransf orraation raethods require to convert the 
general eigenvalue problem: 

=: I^M into a special forra 

! aJ X = X X 

and then Involve a series of similarity transformations 
which aperate on the matrix \ a J . This requires the 
storage of the full n x n matrix 1[_A jin the computer. If 
n is large, it occupies sufficient computer space* I'urther- 
more, all transformation methods like divens, Householder, 
Jacobi etc,, have been found to be efficient for the 
complete eigensolution i,e,, determination of all the 
n eigenvalues and associated eigenvectors of the set of 
equations Sqn, ( 3) , presuming that there exist n distinct 
eigenvalues. 
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However, a reasonably accurate dynamic response 
analysis oX structures can be performed with the knowledge 
of only the first fevi eigenvalues and eigenvectors from 
the lov/est end of the spectrum and hence the complete 
eigensolution is not required, Therefore, transformation 
methods are not very efficient for use in deteruilning the 
dynajTiic response of large sized structures. 

Iterative methods like Power Method, also require 
the initial conversion of general eigenvalue problem to a 
parbicular one and thus suffer with the same disadvantage 
as the transformation methods so far as computer space is 
concerned, furthermore, these methods are not very 
stable a.nd require periodic cleaning inorder to obtain 
eigen\ralues higher than the fundamental one. 

The iterative method described and used in the 
present v/ork has the advantage that it does not require 
any initial conversion of general eigenvalue problem to a 
particular form, ivloreover, the method does not require 
the assembly of (n z n) structural stiffness and mass 
matrices. The entire algoritlim uses the stiffness and 
mass matrices of the elements which compose the structure 
and their order is very small compared to n.. 

Lasbly the method has been observed to be 
numerically stable and efficient for the determination 
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Of partial eigensolution i.e«, the first fev/ natural 
freq.uenciGs and mode shapes of the vibrating structure* 

Oliapter 2 describes the formulation of the 
method for eigensolution, Chapter 3 describes the 
algorith'n used to obtain the eigonsolution. Chapter 4 
gives the results for the class of structures, to which 
the method has been applied. 


« 


/ 



- II 


rOHi/tULATIOR OS' THii' aUTHOD vF SIGI^vlSOLUTlOl^]* 

All iterative method using the v/ell 1 -vTIoivti property 
of RayleigJi quotient is presented vdiich can he applied 
directly to the eigenp rohlem: 

! K j X = ^ 1 M 1 X 

It uses the property of the Rayleigh quotient: 

-y-*- r 

R (X) = ^ X j that it equals the 

X j i<I : X 

eigenv:ilue v/hen the eigenvector is substituted into it 
and that it is stationary in the neighbourhood of an 
eigenvector, 'the R-ayleigh quotient is minimized to obtain 
the lowest eigenvalue and the associated eigenvector. The 
minimization is done numerically using the conjugate 
gradient method. The intermediate eigenvalues and 


* 


The method was developed and is described in detail in 
reference no. 3? or completeness of presentation, .it is 


briefly described here. 
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the aerjocisi'bed eigenvectors are ohtaiiieci by using a 
gradient projection scheme for constraining the minimis- 
ation 2 0 \rc h to t he sub sp d , c e - o rt ho g onal to t ho p r evi ou sly 
determine cl ' eigenvectors , 


'fhe advantage of tho fo omul at ion 



X ( M ! X 


is that both tho ivinorator and the denominator^ as v/ell 

as all of the other q.u antities needed for iteration 

procedure for all of the eigenvalue desired, can be 

computed v/ithout having the assembled ' K • and ! ivl \ 

matrices at hand. Since, the immerotion is tv/ ice the 

— > 

strain energy for a give X (the displs.eement vector), 
and tho denominator is twice the maxm'om kinetic energy 
of t he s 'b rue t u r e , the sumin at ion of p o t e nt i al and ki n et i c 
energies of tho individual elements of the discretized 
structural model gives us the numerator and the denominator. 


2,1 f ORLiULibfIOk Of I'Hf f ROuLK,.; 


The eigenproblem can be written as** 
^ \ 

\ [ K ! - ^ i M ' d = 0 
(' ■ - 


If X xs its solution, then b X is also a 
solution for any non -zero value of tho scalar b, and thus 



the eigenvector correcporidins to any eigenvalue is 
arbitrary to bho extent of ..i ycalar multix:)lier, I’he 
R-ayleigh rjiiotieiit ie not defined o.t the origin and thus 
Care has to be t,:.lven v/hile minimising, dhe redundant 
degree, s of freedom v/hich prevents tiie doter^ninati on of 
absolute muignitude of tb"' components of the eigenvector^ 
Can be eliminated by an arbitrary noima],ization and the 
simp 1 G ,st normal i 3 at i on is to 3 et any n on - zero c ovaiD onent 
of bhe eigenvector equal to one, !^he value of the 
layloigh quotient is bounded by the lowest and highest 
eigenvalues of the physical system. 


'thus, bho minimization of the Rayleigh quotient 
will yield the lowest eigenvalue, 'Ihe minimization 
problem to yield the lov/est eigenvalue can thus be 
stated asi 


ai'ind Z 

a (x^) 


such that 
^1 ^ 

1 1 ^ 1 


-^1' 

is minimum, subject to, = 1 where 

is the norraollziiig or roferoiice component and e is 

y. 

a vector with its qth component as one and zero else- 
where (i,e,, is a unit co-ordinate vector for the 
qth co-o3:dinate) , 
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OnCG tile longest ei^env^-luo = A (X^ is 

known, tin next liiglirir or (oecond) eigenvalue a/nd its 
associated eigenvector can bi- do be mine cl by -posing a nev/ 
IB inirni sat i on p r o bl ei:i « 


Find X 


R (X^) 


X2 j such that 



is rn i iii. * uib , s ub 3 ec t 




2 


to 

1 


r'i r . "x 
and X2 i 

in bhe sub space defined by the oJoavc constraints, 
the Rayleigh q.uotiont takes on a unique minimuni at bhe 
eigenvoetor associated with the second loviest eigenvalue* 
The seco'iid constraint represents the imposition of the 
M - orthogonality condition betv/oen and ^2* ^eonietri 
cally speciang, these constraints merely restrict the 
position of vector space in which the seoach for the 
second eigenvector Is carried out and in tliis restricted 
sub space, R has a minimuni corresponding to 


TieterminatiOvi of the third and siibsec|uent eigenvalues 
and the associated eigenvector's is accomplished by solving 
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a sequence of proble''iic siinilar, a.8 for the det emimtion 
of 2nd eiffenvalue. dhe only change Is that each time 
^^fie a^dditional equ action of constraint has to he imposed on 
the mill Lmizat ion proiDlem to satisfy the condition that the 
eigenvector currently helng sought is U - orthogonal to 
all or 'f e previously deterrained eigenvectors* prohleDi 

of determining the 1th eigonvtrlue ( 2 ^ 1 ^ n ) 

Con thus he v/ritten 


i''ind X 


( q ) 


‘Subject to e 


, sue h t hat 


-*‘P r n 

illlh 

..yji 7 "j 


= 1 


_.m - 

and X^ 


I -> 

^ X. 

X 


- 0 , 


i=1,2. , , ,1-1 


vmere i = 1,2, ... ,, 1-1 a^re assumed to be 

Id own when the 1th eigenvector is being sought, 

r . . n '->• 

lienoting; [ ti -j X^ = the above constraint 

equation c?in be wxdtten as 




X-i ^ . 0 , i — 1 , 2 , , 


* • > 


1 - 1 



OHAP — X.lX 


i'.i Hi Xii X ^ Ji XX Oil ~ Jill G OKI 


Xhe gradient niethod 3 arc neualXy more powerful to 
minlrn.ize a function of n variables as t’eey uGe the local 
infOiTaotion about tiie r-ita of change ox function v/ith 
r e sp act to t he c h ting: os In t he v r .1*1 ohl e 3 * Hence one of t h e 
gradient i^iothocs is being used here for minimis, at ion of 
the Hayleigh quotieeit, Xhe -Rayleigh quotient as a 
function of the 11 variables * ^11 ^ 

diff erontic-hle, '.aid its grodient vector 



_2 JC 
' li 1 X 


( X ; h .( X) 

(xta] 2)2 



(' ""t I' "i — > 

2( I .K .__X y,!) > X ) 

( X^ M i X ) 


is e ra s i 1 y c.oin y u t o d . 


Xhe steepest descent method v/hich is often forced 
with the c Oiivai’geuc e difficulties is modified for rapid 
couvergenice by conjugate gradient method, 'Xhe ^llatoher 
Reeves Conjugo;be gradient method has been used here, for 
reasons of its simplicity and its modest storage demands. 
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T‘he algoritliii is as lollov/s: 


3,1 lijiil -jD iiiiiHj \/ -LJtj i'.'iiij j.j iOi) j 

5,0 I'A ’ ‘AK 0. 5a’ A ent^si. 

'i'hG coTivergenCe dif L'icultir-s oi 'uhc steepest 
descent (gradient) laethod caa bo great j.y reduced by a very 
s imij 1 e I'Ti o di- f i c j t i o a v-i fi i c h c onv e rt s j. t tot lie ^ on J u gat e 
d radi ent -method, 

dere, tlie entire algorlt Ixa reduce to 



arl)itrar'/' 


r; 

'O 






0 


0 


0 



X. + 

1 




i+ 1 




§. 

X 


Clearly from 
1 ine ar c om inat t on 


fchis definition of 
of 3^^.^ and Sq, S^, , 


it is a 
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and hence j it: is a liiiefir cornbinabion of G-q, G-^, i.*, 



^he conjugate gradient I'lethod is a rt,wpidly convergent 
techni(.[ue, suitable for use v/hen the gradient of the 
function is readily computed. It generally far surpasses 
the steepest descent met bod excel) t in unusual cases and with 
Careful applicationj it is one of the most effective 
minimizuotion techniques. In px'Ogramming bhe method, 
attention is to be paid to the questions of the accuracy 

•¥r 

required in the determination of a, a strategy for 
restarting and a proper scaling. 

there are more pov/erful methods, but they are 
generally impractical in large complicated problems. 

the Conjugate gradient method is the most 
efficient minimization technique as the problem size is 
increased, iuspite of its wealcer stability, 

■I'he problem of luiniraizing the Rayleigh quotient 
function to find the lowest eigenv^iluc is similar to an 
unconstrained minimization problem and the algorith can 
be directly applied^ the use of conjugate gradient method 
for finding the iutennediixte eigenvalue is possible only 
when the minimization of the Rayleigh quotient is 
restricted into a subspaCe of X in vi/hich the constraints 
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I'l e. = 1 .-md X '- MjX = o, i =1,2, -1-1 

is continuously satisfied. I’lius, in order to insure 
that tliG Fjearch is carried out in the desired suhspace of 
X, it is necessary ( l) to start the iterration with a po- 
int ii'i that sub space and (2) to project the gradient 
vector g, onto' tliat subspace. buch a projection raatrix 
i s g G n o r o.t o d b ol o\v ■ 


3.^^ PROJEOflOlj LuiTilLX: 


1 


J-ietl P bo a matrix v/bich ha.s the property that 
for any vector ’//, the vector 


l T — >■ 

\i = » P i w 
p - J 


Satisfies, W 2, 

P j. 


= 0, i = 1, 2, , . . , q 


v/here S-^i= 1,25 . , , ,q are q linearly independent 


vectors , 


■the abo^^e -bcjn. can also be written in matrix 


form as 


r - ^ -> 

! i'l W =0 


(qxn) (nxl) 
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where 




♦ 


.'nxq 


1 



z 1 


In otharv/ords , the raatrix operator 1 eliminates 

^rom th: vector W the non~orthogonal cojjip oaents, thus 

— ^ 

giving the vector W which is orthogonal to the suhspace 

— h 

Spanned by the vectors i = 1,2, , ♦ , ^ g, I’his 

idea is differently expx'essed as: 


\1 

P 


= V/ 


1 

2 Z . , or in matrix form as 

i=1 ^ 


V/ 

p 


I . _ \ 

W - L ^ J’ U 

( nxq^) ( qxl) 


where the components of vector U are Uj^, i = 1,2, , , * , q, 

r 1 1 

P re-multiplying above Eqn, by W J , we obtain .from 


P 


= 0 


f = [hJ® i - [n] ) u = 0 


u ■ = ( [tl j * \ h] ) '' [ ffj' 


w 


bus, we obtain 


% 


V 


= w -[^wj M)"'' ?/ 
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(K W) v/ill be a ( x q) symirietric matrix, and is 
norLSinqrilar since 1 i'* 1 is a ( u x q) matrix composed of 
q linearly indexDendant colui-ins. 

HGiice, its inverse exists, 

The projection mabrix^ by comparison is obtained as 

W = ^[1] -bO ( )"^ 1^0 "'j 

S' or determining the 1th ei gens olut ion, a projection 

matrix 

M i WM"i] >-'K3'] 

(nxl) (ixl) ( Ixn ) 

where ^ e.^ > "^2^ * * ‘ ’ ^l-l J 

will project the gradient vector g, on to the subspace 
of constraints already defined* 'the colujnn vectors of 
the matrix linearly independent. i?has 



One thing that warrants attention here is that , 
since we use gp instead of g , will the conjugate 
gradient method actually converge? The answer is 
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afiriiTaativv? and it can be Ghovm by constructing the T^ag- 
range iVuiction 1‘or the un const rained miriL-jii nation problem, 

Siria? BY STEP PROCEDURE: 

5*5 (a) Choice of a Starting B oint; 

A the iterative methods should have a good 
s t art 1 ng p o i nt , ot he rwi s e time is u n~nec ess aril y v/as t ed 
inside the minimisation procedure, A starting vector as 
a unit vector (the jth element of is lying 

along the co-ordino^tc axes, gives ;an awfully distorted 
mode shape in thedynamics problem, A set of random vectors 
proved to be superior contrary to the expectation and in 
this v/ork they v/ere used as starting vectors, i'he 
startiiig point for the search of the lov/est eigenvalue 
needs to satisfy only one constraint = 1 

w Viic h is 1 3:i V ial 1 y s at i s f i e d by di v i di ng t hr ough Out by t he 
qt h c orp on en t . T he m i n im i z a b i on al g o ri t Imi u s i ng 
conjugate gradient method generate a sequence of vector 
which, in bhe limit, tend directionally to the minimum 
eigenvalue onthe search space* Ihus the qth component 
of the vectors so generated ha-Ae to be maintained as 
unity, so jis to satisfy the constraint Aqn, ®q~^ 

c ontinuously in the space, This is achieved by setting 
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tile cit’a oonpoiieiit of tlio gradient vect-jr at tlie particular 
point ecu.:'! to zero. In other v;ords, no 'move' is made in 
the I'ch diroctioii oj. the search space. I' he starting point 


for the searcli of second eigenvalue hcs to satisfy o.n 

J? * 

additional constraint hq_n. this is 

achirured by 8chuidt orttiogonalisation, 

r( o) 

let -h he some iriitiod, estimate, 

fheroCore, ^ ^ 

T 

satisfies the aoovs constraint if V 1, 


Since V” ^ = (M \ where is the eigenvector 

corresponding to thf; fii'st (lov/est) eigenvalue. 

fhe starting p oirit for the search of the subsequent 
eigenvalues, is obtained by passing the initial estimate 
through ru projection matrix, where the projection for 


, the 



( nxl - 1 ) 


folloV'/lng 




.1 ( b) I'Unci'IUn ^/.Jjii^,fl01: 


It is necessary to have an efficient routine for 
function evaluation, in order to avoid the time YJhich 
would be otherv'/ise wasted inside the minimization procedure^ 
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(^0 


Can l:e v/ritten 


K (.C) = 




r 

Z 

i ^1 


“r 
L ^ 


L 1 i 1 


L^-iJ 


V7 lier e r i 3 ■the nuiti her of d e s c r e t e el evii c n t a j i an d I' ^ j_ \ 

are, respectively the stlfi'aces and >iiass matrices of the 
itli element of the discretized structure aad is the 
disi:slaeei'icnt vector of the ith el era e at correspoadiag to 
the genersiiized displacement vector 


Ihe size of element stiffness and mass matrices 
are reJ.utively much smaJ.ier tiian the size of the assembled 
stiffness an d mass m at r i c e s o f 1 ar g e c omp lex structure* 
l“he stiffness and mass matricas of various Id-iids of one 
di e'M e isi onol. st me tur al el eixie nt s are gi c a in Ap p end i x - , 


3*5 (g) dp-adient Avedyati,^ 

'the minii'7isatioi'i algorithm require the evaluation 
0 f gradi cut vector at , e ac h c yc 1 e of t h e i t e rat i on ; 1 he 

gr adi e nt vector, g , of t h e e i g h qn o t i ent f unc t i on is given 

earlier, ^'or determining the lowest eigensolution, the 
coiTiu onoiit of the gradient vector corresponding to the 
norraadizing o opponent of the eigenvector- is 


set equal 
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to serOy I'l order to coiicinn.ouoly satisfy tho constraint 
impooed die to in.e norraalisat ion oP the eigenvector. I'' or 
t he i lit e ra e cii at e ei g e ns ol u. t i o n , t lie ’ p ro j e c t i on mat r i x ^ • 

ic uGed to pr ject this gradient vector, g, on to the 
proper sue space of search, ii - orthogonal to tlie 
viouslsT determined eigenvector, Ihs component of the 
gradi out vector, gp , c orre sp ondi ng to t Ue no rrnal i si ng 
c ornp o n ent o ;r t he ei genv e c t o r turns out to be s e r o au t om at i c a- 
lly, but 0 . small number generally ao pears because oi 
roundoff errors and this is removed by setting: that 
Gomponont equal to zero. 


• 6 ( d) .^yal.iicfbioji ^of. .S tep.. L^engt]!*. 

l^iice a direction of uovo has been decided, we 

mu s t de t orm 1 ne a sot ha c t n o . ^ uiio t i on i s minimi zed in that 

direction* d'h.i;i 3 iii' problem of finding the step length 
is essentially tho linear search * 

/problem' "vhicli requires the determination of the s-long 

3^ tljrou-.<h at v/hich uhe value of the >tayleigh quotient 

function 


■I t ( X . + Ci . S .. ) s= 
1 i ;l 


/ W -♦ N 'I' ^ ^ 

( X . + a. S . ) 

V T •! n ^ 


. 1C ;(X, a, 3.) 
11' u J ' 1 11 

(X. + a. S ) ■^' r li 1 (X. + a. 3 ) 

1 ri LJ i li 


is a :'iinimuni, i.e, , 
dH/da^ 


a. =3 
1 1 


= R ) = 0 

1 
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In the general arovyic:., no exp re ee. Ion is available 

■V'" 

t o dot G na i ne , so a .'l i nt e xn ol at i on ap o non. c li 1 3 adop t e d , 
iiowevex'j in then articular prohlem 01 tli:- la/leigh guotirent, 
an explicit expression in cnn be generated, which has 
the fonn 


nL(. . T* va . t = 0 
" X X 

"'•“I r -1 T-i' r ■ 1 \ 

i^vhere u = L ^ 0 ^ ^-^1 I J 

- (x^’ [ ids.) ( shin s.) 

^ X ~ ^ X ^ 1 s -i 


= (s'l [k] sp (x'dgx.) - (xbicl s^) (s'J; [m] s.) 


= (a h'd x^) ~ (x^ Lx] 


I he tv/ o 1 • o* ts c 0 r r e :-3 p ' nd to the m nx Ini -al an d i.ii nim al p 0 i nt s 

of the h.ayleip'h quoti nb in th, direction tlirough 
'1 

■xhe .ainijiiE’l function value; oorrosnonds to a. . 

1 



GHiiPTjliR - IV 


Oi] — Ai'U.) ~ idJiD J-UlS 

In aJ.l the exarnpleri p resented in this ^iiapter, 
the stiffness and mass matrices are talcen from Appendix - II 
Ihe dist.rihuted maso of the sj^stam v/as ij.sed to ev^.luate 
the components of bhe mass matrix (the 'consistent' mass 
matrix) . 7044 digital oo'.mputer v/as used to obtain 

the numerical results and random vectors v/ere taken as 
starting points for the i:iinimiKation algorithm to obtain 
the eigensolution. 

4 . 1 .'fiXAIvil’liA ( 1 ) : 

As a simple application a slender solid cantilever 
rod shoeni in I'lg, ^ l) is considered* I he first tliree 
eigenvalues and eigeU'/ecto.rs of this simple structure 
^v;ore obtained by the isaccessi/e minirnir.ation of the iiayleigh 
quotient in the appropriate subs^ace. 

I he Cantilever rod shOAvn in I'ig*( -j) is divided 
■ into three elements and each of them is modelled as a 
genersl planar bean element having six degrees of freedom. 

I he cantilever Iig.( 1) thus has a total of nine degrees 
of freedom, Ihe first three eigenvalues and thGir ' 

r 

associated mode shai3es are given in lig, 2(a) to 2 (g), 
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^he -l-^able ( 1 ). f conipareG the computed results v/lth actual 
theoretical results obtained from closed fonn analytical 
exp re s si on , 


4.2 Another illustrative exaiiiple considered here is 
the portal frome shov/n in I'ig, 5, Aa-«h meuLber was 
model-led with a single beam element, I'hus the structure 
has four degrees of freedom. A complete eigensolution 

was obtained by successive minimi zatioa of the hayleigh 
quotient, Ihe time taken was 24 seconds. 

The computed frequencies are given in Table 2, 

4.3 i^inally a skew bridge structure with 12 degrees of 

freedom shown in ^'ig, 5? has been analyzed to obtain the 
first four eigenvalues and the associated eigenvectors 
by successive minimization of Ihe analysis is 

Carried out for a grillage v/ith two longitudinals, fixed 
at the extreme ends and their cross girders are assumed 
to be simply supported over them. The torsional rigidity 
of the members is neglected. 

The computed velues of frequencies obtained are 
Compared with calculated values of the frequencies and 
the e:q)erimental values of frequencies as given in-def. ( 1 1 ) , 
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and arn indicated in Tadle 3. 

’-^lie first mode shape is plotted and is sho\m in 

5'ig. 6, 

file total time taken to compute the first four 
eigenvaJ-ues and the arjsociated eigenvectors is four 
a'.'' I > r o X i ra a t el y . 


m inut e s 


TAiJJjiJ 1 


It'requenc; 

r — 

/■ ciiai Values iii 

'compute (I values ir 
c;s/sec . 

Ao. 

Co , / GfiC . 

1 

3 . 1 c 

5.099 

2 

51-599 

31.799 

3 

f — . 

CO 

♦ 

CO 

CO 

S8.99 


2 

^;0AiPUa'E13 PRtS;QUiiii^iCIjiAS I A OS. /SAG, i'‘OK POAl'AL 

Akiil; 

1 . 42.99 

2 . 241.23 

3 . 833.23 

4. 1559.72 



'-o-’ii, iiL jii 
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I’ renneacy 
I'-'O. 

3jL c 1. f 1 1] 0 d 

Value of f. 
iu Gis/sec, 

— 

^7.[j ex'lnien \'>ril 
volue of f 

in cr:/see 



Computed 

value 

in cs/sec. 

1 

11.136 

11,222 

1 

11.605 

2 

12.72 

12.53 

12.80 

5 

27.82 

27.98 ■ 

28. 30 

-w 

4 

3: . 44 

34.21 

36.62 


i. . . 







Sjjiili'iDiiiR 


Cajj! 'XXXjjiv ER RuIj ( 9D , vJ ,i^' , ) 

L = 2^5,0 in. 

3 = 1,0 in. 

6 , / 2 

^ = 00,0 X 10 J. us./in 

y _ 


0.2G Ibs/iii^, 




Sect j o m 

S £ c T 1 O w 

Ax- A c-c 

b- 


PIG-. 3 

PORTAL inLAiL j 

POUR UPGRERti OP PREELOA StSIIiu, 

L = 

25.4 cm (10 in.ohes) 

E ~ 

2.1 X 10° kga/cni^ ( 30 .O.S 10*^ llis/ln^) 


7.81 gm/om^ (0.2a Ibs/inch^) 








51 








APi-'jiljDIX - I 


JjISx of Ri;;FJi;itFi'[CES 

1. Wilkinson, J.H., 

( 0:^::C ord , ^1 ar endon P r es s , 1 9 G 5 ) . 

2 , Ant Uony dal s t on , ^ s e^. In N umerical .IncvlXSA j 
( iivi cO raw H il 1 *9 o . ) . 

3 , R , E . P , Bis hop , y.1Jbrg^tXQJi ri dg e U ni v e rs i t y Press), 

4. Pox, R,B , and Papoor, i'vi.P., d ..Minimization i-iethod for the 

S qlu^tion .of ,t.he .^.igej^.i’pbl.eni A^.iP.irig ,i?V_S^tru^tarcvl^ ®j;n£gr]i.i,c^s_j. 
aPePL - PR - 68 - 150, Paper presented at the V/right 
Patterson Airforce Base 'Conference 1968. 

5 . J .B , P rzemleniecki , jJ^Ueory .Qf B t rue t u r al, Ap al y . 3 1 s , 

( McOraw-Hili Go . ) , 19 68 . 

6 . Har ol d , 0 , M ar t i n . , Pnt r ojf u^c^t i on ^t^q s^tje^i i'£§ .U^Qjds^.i? J 
St.rLiAtdi3;.a9-^..Ana;Ly^?IX^^^^ 

7 . Rubinstein, M os [ j e , P . , M.atrt;^ G_omp.uter,.,An,gJo(sis.^ 

, ( IT ew J e rs e y ”• P r en t ic e -B all , 1966). 

8. Gere, J/M, and Weaver, Willi am, Jr., Analysis of Pramed 
Structures (Rev/ Jersey? P. Van Rostrand Company, I 965 ). 

9. Jacobsen, Bydik, S,, and Ayre, Robert, 3,, Engineering 
lxl>xa,tXom .v{ivt.R 

(Rew Pork: PicGraw Hm Book Go, 1958). 


34 


10, Hu r t y , W , 0 . , and Hub i n3 1 e i n ivi o s he , • , Dynanic s._. q ^ 

St ru c _t u r es , ( i'l e > v -Jersey" J' r 1 n t i c e al 1 , 1964). 

11. A rm strong, I . D . , .4' l ,..e _ i iatu ral J'requencl e s . of .'J;r i 1 1 a^e s , 
(international Journal of i-'iecl anical Sciences, 1968, 

V ol , 10, pp ♦ 43-55). 


12, Armstrong, I.D., .^aiUQS^qf Oynap^io Stifxness AAA. 


y (Sfvil J^ngineex'ing J^epartment , Heriot-Watt 
University 1967) . 

13. -9,J, Allman Eind U,M. Brotton, 'Systematic Relaxation 

iletUod for the Uetervaination of tine formal Uodes and the 
Uaturol frequencies of Vibra.ting Systems' , (International 
Journal, of Mechanical Sciences, 1968, pp. 129-141). 


14r Mahesh, P. Kapoor, Automated Optimum Besign of Structures 
Under dynamic Response Restrictions, (A Thesis submitted 
to Case Uestern Reserve University in Partial fulfillment 
of the requirements for the degree of Doctor of Philosphy, 
1963) , 

15, Pox, R.D,, Qotiniization Methods for En gineering De sign., 
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fi/lASo iviAX'iLLOijrj ii'Oil VAfaOUS 

Ojj 6 X iriv J 0 TU Iu-lXi XlLi jA. l jl;lt XS 



A Lffan iiljAivIiiirX GOllSID^RIi'^a LATHlKAL 
VIJ.RAXIOKS 



A Grittr^ JIJH jlL k^X^AOl’j 3? iliTi mVj r^ij'J'T 


(12 D. O. P) 

Note* Double Arrows tno+cate Rotations 

AnJ3 the SIN'a.uc AR,ROW INJ)ICA7ES 

D i S P L/kce mints. 













STIFPM5SS MATRIX FOR GEI®RAl BEAM ELE?.IENT AWIi FOR GRID 
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- III 


00j)M?UTj5R PROu-fOiM 

his appendix conialna a description of each of 
the i'ortran subprograms which makeup the Oomplete program 
for the free vibration analysis of structures .composed 
of elements described in appendix - II and the listing 
of each suloprogram is appended £it the end. 

Jir .2 f rqgxsyplL 

Ihe complete program for the free vibration 
analysis of structures consists of (a) a main program 
v/hich reads the input data, evaluates the lengths of the 
various members of the structure, generates the starting 
random vector a.nd prints the various input data. and 
( f ) the other subroutines used ere described below* 

^ ( a) of ^iT,'brouti ne» KMK 

( b) Aj!§iWAcnt LD?, RDF, .RI)F3 ^ KL, ivED, KR, 

HR, 


ivl : Total number of elements 

IjI)P : local degrees of freedom of an element 

RDI' ; Number of degrees -of -freedom of an element in 
System co-ordinates. 



NDJi'S : i'! amber oj? degrees of Ireedom of system 


ICj ; ii^lemeut stiffaess matrix in local -'Coordinate syst 

ML ; L’lement mass matrix in local -c oordinate system 
idi : Ele':^ent stiffness matrix in system co-ordinates 
mil : Element mass matrix in system co-ordinates 
T 11 : '.0 j;an 3 f 0 rm at i o a m a t si ic f o r an el em e n t 


(c) C_o;‘.ir'\ort Va^r fable sd. 


L, L-IIX, lilY, MI2, 


RO, G, ALPHi, 


GXj CY, CZj Pf], Q3^ IIPROL. 


'• Modulus of Elasticity of the member 
1 • Length of the member 

MIX ; itioment of inertia of the member about x-axis 

MIY ; Moment of inertia of the member about y-axis 

MIL : iioment of inertia of the member about z-axis 

RO : Weight per unit volume of the material of the 
member 

G- .* Modulus of rigidity of the memher 


AJjPILI : Angle of inclination of the member 


GXj GY and GZ are respectively the L and Z co-ordinates 
of the various nodes of the system respectively. 

PE and QE : are respectively the left and right hand 
ends of a member. 



'M 


rip ilOB : X rj a counter in di c a in ^ t he t yp e of t he p rohl em . 

indicates 

£ that the system consists of truss 
olemonts, 2 indicates that the system consists 
of either the general planar be^oit element or the 
grid eleiiiont for inyjlaue vibrations, 3 indicate® 

thot the system consists oi grid elCMnents for 
leteral vibrations and 4 indicates that the 
system consists of space frame elements, 

2 ( a) AnpTO ,.qf SpbroptihO' ZA’U-il 

( b) Argument List: x'.'i j HbP ^ SR^ K.Rj 

'iost of them have been defined earlier. 

SR : Assembly control table. 

( c ) SBIEf' ^ aA j J3B ^ NUldER, DJSR Oii'I ^ VRG ^ 

BX, ivillA iViXRO, I'fER, ORBdA, VALU, 
WBC. 

SDIRP : Scalar; denotes convergence of Rayleigh 
quotient function; (rAAId) 

AA : Vector; denotes where ^ is master 

stiffness matrix and X is generalized displacement 
vector. 

BB : Vector; B enot e s | v/here 1 is .. t x: 

mass matrix and X is generalised displacement 
vector. 

FUMER : Scalar; Denotes 
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jOEWOi.i ; Scalar; l^enotcs 

VEO ; ^oiAiiter; l^euotec the number of eigenvalues— 

eigenvectors being -searchecl. 

BX : Matrix ; the ith column of BX represents the 

vectorj^hi J v/herej^ h ^ is master mass matrix 
— > 

and rex^resents bhe ith eigenvector, 

HI I'T : Counter; Denotes t he c u it ent c ornp o n e ri t nurn her of 
the eigenvector v;hich is set to one. 

MlhO : Counter; Denotes the old component number of 
the eigenvector v/hich was set to one, 

ITER : Counter; Denotes the number of iteration of the 
Conjugate gradient method. 

OH EG-ii : e c t o I' ; t he c omp o n e n t s of 0. EG A r ep r e s e n t the 

frequencies (square root of eigenvalues). 

V^ILU : Vector; the components of VALU represent the 
eigenvalues. 

VmCO ; Matrix, the rows of VmOO represent the eigenvectors 
NVEC : ilumber of eigenvectors used for modal 
t ransformation. 
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( d) 'i"o compute the eigenvalues and associated 

eifpnvoctors* 

3. (a) pX J' jU'iDG 

( t) ) A errt ...L ip t f or t hi s su b r-ou t i n o t he ar gum ent list 

has already been defined. 

Purpose i lo compute the function (Rayleigh quotient) 
value o.nd to evaluate the gradient vector to the Rayleigh 
quotient zit any point, 

4 ( a) i'^oino of Subroutine* ^ROD 

(b) / .1 ® defined, 

„?n. 4 TP>oseJ. ‘^o compute the product 

vzbere 

and '.V is any vector in vsysbem co-ordinates, 

5(a) R smi e of S u br ou t i 

9 

(b) Ai*^ument J^istj Already defined, 

? Lirp o s e ^ o obtain ( [. R J ^ I J ) ^ 


from 
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6 . (a) ^ of ^^^ubr 9 ut • lOUil. 

( b) *. WD5'o , XRy LBi3 

NBES : li umber of degrees of freedom of an element 
in system co-ordinates, 

SR : Assenbly control tixble 
X : Vector in system co-ordinates 

KR : Bleiiient stiffness matz’ix in system co-ox’dinat es , 
idR ; Element mass matrix in systeoi co-ordinates, 

LBf : Jjogical; Ebb = ,!rRUE. provide the normsdized 
vector bB = (Euclidean norm of the 

vector is equal to one) , 

(c) 0 OiiUa on y ar i abl es.l AA, Y, 

AA : ector; Benotes where ^ "j is master 

stiffness matrix and X is generalized displace- 
ment vector, 

BB : Vector; Denotes f bJ where is master ' 

mass matrix and X is generalized displacement vector, 

X i Matrix; the ith row of the matrix denotes the 

divsplacement vector of the ith element corresponding 
to any generalized di placement vector X in system 


eo-ordinates 
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^ V ° ° ^ ^ ® ^ ve ct ors ^ ^ auid iVi 7^. 
where ^ Iv'jand [.,!] are respectively the master stiffness 
and raass matrices of the stractare and ^ is any vector 
in systoiii co-ordinates, hote that master stiffness and mass 
matrices are not asseobled and the desired vectors are 
obtained hy the use of only eleuienb stiffness and mass 
matric os . 


7. (a) iUv.P, Pf . ?nbrputine; TfiSlS 

( b) ■;Vrf<;uiaent ..^_is^_ Al ready defined. 

( e) '^‘juiaon Yariables^ ^ NDFS, ilVlG ^ X, BX,^BX ^ VliC, 

^ .f VT-fV PPP: 'j^ests for vector space, test for gradient, 
test for movevector. 


listing of each of the seven fortran subprograms 
which malcG up the complete program for the free vibration 
antilysis of various types of structures for IJi»I 7044 
digital computer follows in the subsequent pages. 
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DFTC f/lAIM ' ; 

Integer pe.oEjRdFsSr : 

RPAL L ; 

real MIX9f/IY,MU • 

FORKAT{2X!»10E1K6) ' 

2 Format ( i6F5„ 1 ) . ,0:' 

FORt'iAT (8 ('ix»r:p« ]. ) ) 

FOn;>AT( lXs23[3) 

format (//93X!. JJ-i’Oc OF OiF,lBb;.S = rvr-'?2X5^. . OF ^0'DFS::* , 1 3 ^ *LOCAL .! 

UpOoF = ■J^vI352Xsifr NC^ OF OcGRiliSL 01' Fi?,LLDf'‘ =■!.' s 13 s / 0 s ? X > NOo OF M 2. 

•ttcIGcNVALUtS And EKjFMVLCr:'S R£:OalOLO 9 13 ^' 2 ^ 9"^ ..pfOu - ^i,l3//) . G,' 

Format ( iHo^ 2 X 92 7 H asslmoly Cl/i.'Tkol .'atrix o^oz/^zx^n.: i.'c^uEp. hou, -i 

l^x 9 19H DEGREES OF FREEEX/i ? / 9 20 ( 6 X 9 I 2 9 1 r’Z 96 t U 9 2X ) 5 / ) ) i 

Format ( 1 H 0 <;, 2 X»’ 27 H ASSEFibl.Y COt.TROL I’ATRIX SR 9 //J 2 X 9 IIH DE-LXR NO., 9 ,'. I 

0f>X>19H DEGREES OF FRELDC.; 9 / 9 20 C 6 X 9 1 2 9 1 ( 9 A(I^92X)9Z)) E • i 

, for'mt ( 1 HuvZXpi'TH assembly Control .,ATf:ix sf; 9 // 92 X 9 I 1 H oEmucr N'j<*9 ) : 

^t5X»19H DEGREF^S OF FRELT00ri9/9T30(6X9l29TrXrl2(l292X)9/M 
■6' FORYl/i't{,13 aXisFS.l ) > 5 

9 F0RHAT( 26F3'u 1) ^ j 

!3' FQPMATt 16F5-1 ) ^ 

data PT/3w]^12/ : 

Coy ,'D^7Gf:AjM/\ZFC0 j 

COMVo'i/KAUL/ 

CoH'"OM/! !AR VIFF j 

Cot-riOr- /SINGH/M s»RDK 9 LDF 

CoKZ'’0;‘VFMT I i.AK / fWF.'C ii 

C.0;}:-.0.'- /level/ c 9 A 9 L 9 . .IX^mI Ys.i iZ 9i';09Gv *.i.PHA 9 CX 9 CY 9 L 7 . ,PL 9 ('t 9 1 'iPROO • 

.Dr<VE/'l$rOfF F( -/()) sRO( So) ,CX( 30) 9CY( 50 ) C.c/!'^’}) 9 FF. < 50) I 30) 9 6 '‘< < 30 » 1?F.aUP 02 OO 
*) sAt 50) >0(50 ) 9VEC0( 10 , 100 ) 91 .1 X ( 3.; ) 9 ! < I Y ( b ) 9 . 1 1 / ( 30 ) ?F( 50 ) . ALPHA ( 30 ) KaUO029(S 
IFF= 0 i 

- INPUT data ^ 

R E AD { -i 9 3 ) ' 9 w 9 L D F 9 liD F ^ NUF S 9 M VE'C ■> MR RL^£ ■ 

READ r'ODAL CO-ORDlMATcS , =; 

ReAD(5,9a) (CX ( I ) 9 I^LsP) , ) 

READ(5998) (CY( I ) 9 I = UN) J 

READ Cj 9 90) ( CZ I I ) 9 1 =1 rN ) ^ 

C READ Sl-CTION CHArBCT ER 1 ST I CS i 

REAIKS,^^) ( A ( I ) 9 1 = 1 sM) .[<;AUv.*035! 

KE AD ( 5 9 2 ) ( Alpha I I ) 9 I = X * l*. ) ; 

REAPCJpZ) (MIX ( I ) 9 I-: 9rU ' 

READi^iO) ('‘lY(I)»I = lj:r) E 

• REA:).( 372 ) (HU( I ) I-T r’H ^ 

C READ MATFRI/Yl P'-'CREfn I E3 . . , } 

■r€ADL 5 *S) : (,E.(rDa=l.H) • • 

rfad( 5'«.»3) ;(G(’ r)vi-i>K). , ■' .•.•5 

. RE:AD(>!r3) ffe6(i)a=TaY) 



RE.-.D 

C(Yl MC 

IUEMCL ATIdX 

READ( 

S d> ) 

(PE( I ) d = l 9 ' 1 ) 

re;'m 

Id S ) 

{ 0L( I ) 5 1 = 1 9 ^ 0 

Rh AO 

ASSd’ 

LLY Cd^TPCl. 

RC:U{ 

1; dd 

( (SrMdj) 5.1 = 1 


CC r.C!-tO PPIMT -- — 

wra T t£ ( 6 > 6 ) . p • . 1 9 L j. r i'.‘ > :M V ..C - . ■ 'i-‘ 

WRITLtf-)?!) (CX( I ) yCY ( I ) .C/( I) 9 1 = 1 . J 
'.JRITL(691) (-( I ) = 

v/RITt(6vl) I Y( I ) 1/ ( I ) ^ 1 = 1 >1.) 

,v.'!' I TE ( 6 9 1) ( E ( I ) 9 no ( 1 ) ? 1 = 1 .:<■ ) 

WRITE(6vl) (0(1)91=150) 

WR I TE ( 6 9 !:) ), ( PE ( I ) 5 0 E ( n V I - 1 ) 

V'JR I T E ( 6 9 1 -^1) ( I 5 ( 5P ( i j J ) 5 J= - 1 Rl F ) , I = 1 ■> ■ ) 

CC ■' GFO'-ERaTF THE ST^'PTtRC V^CTOi^S 

u6^ 1 ^ 1 = 1 5.HVEC ■ ' ' ’ ' 

■ ' DO 12 'J = i ;HOhS 
, 2 ■ ' VE C0'( I , J ) = ' ; , P,-!DY 1 ( . u 1 

‘D(^’.f)'3 I-l 5 r/ ■ ' 

Rod) = RO(t )/9filoO _ 

; calculate lCmGTH of Amu UIRuCTIO''! CO^jItltb 

JCO P0!3AI = 1. 9i'. 

JT=P7( r ) 

'A ' M ^ 'sGRT ( (CX(RT )-CX( JT) ) ^^*^2+ < CY ( X T ) -CY ( JT ) ) ’'^^^2+ ,CZ ( KT )-CZ( JT> 

■ W --' I T E ( A 5 3 1 6 } ( L ( I ) d = 1 9 .^ ) 

.Call f r • i m ( ■ ■ d df .. vor s 9 s ) 

STOP 

EMD 


TC Kr-;R 

suiir>(^u r IHE x/'rc^p? ’ d) 

1NTF‘-EP Pr’";lddG!-. 

Real u ? l .l 9'-' 1 •; ■ '•i ; l ? ? 2 ‘>r'U- 9 l-o 9. 'O 5i\3 9 ■ 0 1 ^ . .02 « kR 5 yiR » 5 rkL 
■real 'iXv-dY^HlZ 

COY -O': /SIRCK/^- yRL FdLUF ... . m. 

Cor-''0' /LFVEL/ L .;A 9 Lj'‘ 1X9 il Y ‘j’i'll? ii'iOo 0 9 />i.Pl''-A 9 Cx jCY 9C2 

loi^ r. { Lo ) , A , J O' ) , u 1 5 c 1 .. a' y, ( ■; 0 ) , . J Y ( ■;! -3 ) , !• . I Z ( S 0 ) , ICO ( 5 0 ) , Q ( 50 ) , 

•^.AlPHA( s. ) flCKl'jO ) jCYI fK. ) ,C2 ( ' 0 ) ( 12 d 2 ) v'iL( 1 2 91 2 1 9 AR ( 1 2 d 2 ) 1 

12 f-]?.') di-d/ vl 2 ) 9T0 ( d.-.d2 ) 9PC(<3C ) 90C( SO) ■ 

• f OP'-^At ( 1 X s 6E1C •• P y _ ^ P . 

CAiL'CUEATU. ;iL L' 'E ' 'T ' S.T 1 F F SS' MASo .ATRlCfcS, 

DO CCF- ■ ■ , : ■ 

'■ ■ D;o- 76 'i-OF L 

• KC.(K9U')'^ • OdQ. ■ ■ ■ , . , 

H[L.(K.»J)'= ■' , " ■> 

■ =KL''(1.>'2)!s' I ) ^fAH.) / t ( 1 1 


54 


i<L(Ji,2)= !<L(]. ,1) 

GG TO ( 10 ,20 9 30 P^OC 

K L ( 3 , 3 ) = { 1 2 c ( I ) 1 1-Z ( I ) ) / ( L ( I) X ^r3 ) 

KL(3,/-r)= -<1(3 93) 

KL(^Fj4>= XL (3 ,3) 

XL (3911)= C «>3i^-XL ( 3 »3 ) ( I) 


KL( 39 ].;::) = 

XL ( 39 11 ) 

l<L("r,1 1 ) = 

■ -XL ( 3 9 II ) 

KL{^V9l?) = 

XL( ^^9ll ) 

KL(‘3,C)) = 

( 12. O^E ( I Y( I M / (L ( I ) -^^-3 ) 

K L ( 5 , 6 ) = 

-XL ( 696 ) 

XL (.5 99) = 

OoS^KL CL 96 )itL( I ) 

KL(5,lr ) = 

XL (3,9) 

XL ( 6 , 6 )= 

XL (5 ,6 ) 

XL (699)= 

3 u 3 -if I'L L ( 6 9 6 ) L ( I } 

K L ( 6 9 1 u ) = 

XL ( 6 ,9 ) 

KL(797)= 

G ( I )-»f-MIX ( I) /L( I ) 

Xi!(79B) = 

-XL (7,7) ’ 

KMa*«)= ' 

XL ( 7 9 7 ) 

XI (9,9)- 

I ,)i^‘:UY( I )/L{ I } 

ke,(9,,1':-) = 

. ; 'a:^3^h^;l ('9,9) ' . 

rst(iQ,9XC.) 

= '<E (9;9;9|'' 

K L ( -ll » 1 1 ) 

Ly ■( I ) Z ( T ) / L ( I) 

KL ( ll9? 2 ) 

= Col'-tf-KL{ 11,11) 

KL(129l2) 
GO TO 10 

= .XL(l.l9ll) 

KI ( 3 9 3 1 « 

( l2oL^!■E{ I ) “r:-iIZ( 1 ) ) / (L ( I )-i*^^^3) 

KL(3,A)=^ 

-XL( 3,3) 

KL(395)= 

0 V 3-i<'XL ( 3 9 3 ) -?'L ( T ) 

K.L ( 3,6 ) = 

X L ( 3 9 3 ) 

KL(494)= 

X L ( 3 9 3 ) 

KL(A-,3) = 

-XL (3,3) 

!<L(4,6) = 

-KL (3,3) 

KL ( .o,::- ) = 

6o0^s-E( 1 )^tHrZ( I)/L( I ) 

KL(3,6)= 

0 r,3-<(XL ( r-> ->3 ) 

K L ( 6 , 6 ) = 
GO TO n 

XL(i^ ,3 ) 

KL{1,1)= 

G ( I )^:vUX( 1 ) /L( 1 ) 

KL ( 1 9 2 1 = 

-XL(l9l3. 

XL{?,2)= 

X L ( 1 9 1 ) 

iViDI= c 

,( T )i'A( 1 ),/L( I ) 

KL(2,3 )= 


Kl (3,-^0- 


KL(3#b)= 

■ -6oO^XY!;l/L( I ) 


6<,o^O'^nT /L( r )■ 


4oO'Xv'Dl 

Kt(G,b) 

. XI / 'i.h ) 

. / ( L ( I ) ). 






110 


51 

CCC 

81 


90 


91. 



';'|L(356)r: .rij-j^Ld )*lj.0/A^0«;; 

= ,:L(0>3) 

‘‘L ( 4 9 !• ) = -r ;L ( 3 9 6 ) 

«iL(4»C)= -Mdld!) 

I’lL (5,3)= •'.''2^<-13uO/33^s. 

ML (5,6)= ;.02^i9..'' /70 • 0 

-■'iL(6,6)= ;'D2-*d3 

COKTTfur 

DO 51 K=1>LDF 

DO 51 j=K,|..nF 

l-L( J,: )= kL (K,J ) 

ML(J,K)== ‘L(K9j) 

FCi: ./'-TlO,'' UT TT<;kNS--t^ i/’.T d 

DO 81 r=l,LDF 
DO ai J=l,Lr)F 
T R ( K , J ) --= 0 c ■ 


AA= /,LPHA(T) 

JT=PF(I) 

1<T= OP(l) 

Ll= (C^(KT)-C'X( JT) )/L( r ) 

011=^ ,(,CY(KT)-CY( JT) )/L( I) 

Ml - (CZ (KT) - CZ( JT ) ) /L ( I) 

I F ( ABS (,M 1 ) uGT iO r> 9999 ) 60 TO 90 ■ 

D F ; . ’ '■ (SORT ("L 1 2 + 1 2' ) ) 

U2:-::ilr-0,3 AA ) ( AA ) ) /OFMOS'i 

|'■:2 = ’' ' MCOS ( AA) 

N2= (-YHtOl^j-COS (/’■;• ) +1. l*s I ■ ( AA) ) /0i:8i'' 
L 3 - ( L 1 !< '■! 1 I M ( A A ) - M 1 -itCO Sii-.r)) / D E,\1C 
>13= OB^'C' 'f-S 1 M ( pA ) 

N3= AA) +L1’^C06( A/M ) /l.'L' 0; 

GO TO 91 
L ], = n - 0 
N 1 ^ ' ) 0 

12 ~ -'-'’'I '■'C0G( .aa i 

i“i2 = "'90 

N2= SIM(AA) 

L3= AA ) 

M3 = '',*0- 

rOS-lAA) 

TrM30.)=. LI 
TR(,2j.2.)='. 11 

TrC3'>^':)^ . 42 : 

t R ( '4 ? 4 } = ,t‘t2. 

T:r,(,5»3')^' 

■FR(6j6).= , 

GO TO' .(,15i»25¥25^.35) >(M>R0B 

!-U' 

TR.|I'.y5)-« ^'1 


TMX 


57 



TR(2,6)= 

•'!1 


TR(3,1)= 

L2 


T:U3.5) = 



T R ( 5 2 ) = 

L2 


TR (4,6)= 

M2 


TR(5,1)= 

L3 


T R ( 6 , 2 ) = 

L3 


T R ( 5 5 '3 ) = 

L3 


TR(6,4)= 

Vi 3 


DO 92 

K = 1 ,6 


KK = K. + 6 



DO 92 J=l,6 , 


J J ■“ J+6 


9 2 

TR{i<K,JJ] 

1= TR(K9J) 


GO TO 3-'^ 


25 

r (U 1 9 3 ) = 

M 1 


TR.(2,4) = 

Ml 


TR ( 3,1.) ^ 

L2 


T R ( fK, 2 ) = 

L2 


GO' . TO 80 

1 ■ 

T5 


u,i 



' LT' ■ - 


TR(3,.3) = 

LI 


TR(4,4)= 

LI 


TR(l93)= 

J'il 


TR (.2,4)^ 

Ml 


TR(3,1)= 

“Ml ■ 


TR(4,2)=. 

-‘'1 

B.- 

COMTtKtUE 


C 

TRAHSEC^'; 

•I SrlEFMcSS -'.•ATKIX TO 


D064J.= 1 ,1 

LDF 


D064K = ]. ,1 

ROF 


SUM- 0 o'" 



D065K!< = 1>LDF 

65 

SU 1= SI 

L; i-hKL( J ) iMfR ( ,‘.'.KvK) 

6^ 

TS( J,K)= 

SUM 


D067J = 3.,i 

r-‘DF 


D067K = 1,1 

RCF 


SU'V= 3 o 



DO60Ki< = 3. 

,LDF 

63 

SU'v)= 5 

u. +TR (:<.<, J )J^tsi;^K 9 K) 

67 

K‘-v( j = 

5U5- 

C 

. TRA» S-I-Of- 

. ■■,\ 5 S i-rtTRlX Tg .■<LF.'. 

66 

DO 70 J= 

I 9 I DF 


D07"'K = 1 » 

RDF 


3U.K=0t.9- 



D071.K.K = 1 

»LnF 

71 

6 UM"— 6 Ui*'+r- L ( J .9 v'.’s )-^T R t , K- ) 

7'^ •' 

TS( 


.69 • 

• coNnf'u'F. ' 


Cl;GRt> I wATt SYSTEI 


*iCc eCvRJJlMATt SYSTEM 


74 

73 

12 


00 74KK = 1 sLT'F 

SU;''.+ T r J J ) --^T S ( ■< , :< ) 
'^R(JvK)= SU''! 

coi-;ri^'ii 
Riniup^i 
EHD 


‘f^IOFTC FB''-;] . ? 

SUb^:oU*Tl;li: 1 ki' I.-', ( ^ '.'UF 5,--';r C- Sk ) ^ 

^ 1 H Is ’ SUb'F.OUT 1 hli’.iJS Tl IL Lioi_i\; VbCTO'S >-..-^0 CCisT- LSPtJ/:u 1 1 ,(' tlOti'v I 

c IjVKLUcS SY THF i‘tlSATIc::’ V, CF Lh.„ '■ ^YLcilbH \:iUOTIbi'!T Z v, THE ' 

^ ?bO',iEiU^AT'^ '"uTHOb iS USEU FCk i " I F I i’^i I S-^>T u I 

ImTvGLp RbF,SR>VLC ^ 

ivfT.AL : ;U ' E R •> f''Ui i t h- 1 ? I ' U ‘r, a t- UF s ('•P 9 1 sK^ku-'.l^ 

COFFoF/PATHAF/FVEC i: 

'■ .cO'';’-':oii/bblocf./aa,9BEj ■ | 

' Gpw:;?^;; /COBoCi-;/ iIUFFR aDEflOF , 

COi'^F'pM/UBLfCK /O ‘EGA ' ■' [ 

EOi' 0 /jpl:'C:</: 'I-' 9.'i r- '■'9 T T -n ■ 

Cfv O: /D JL^CK/VFC 
C0M'C.>7GBL.0C*'./BX 
Coy- 0'7G.MA ( '-/VECf; 

Co*: 'OF/HARI/IFF 

Oh.Ettsio.'i. SR ( 5u ) sxi 1 ) ab ( :■ ) aS ( ) ax£' ; { loo) vGxE'it icc ) ^aa ( iooi^auoj 

1 ) !,BX (.100 9 1'.J ) a,xr:-x ( 1. Ol y ( ICO ) 9i-- '-P ( lOO ) ( lOr ) ,X''ih ( 10r*'^A(J< 

* » 1 f; )• > P ( I CD ) \/AL C' ( 1 0 ) y Vu. CO ( 1 0 .,100 ) 9 KR (I 2 , 1 2 ) 9 . iR ( 1 2 » 12 ) 9 0) .EGA ( 1 0 0 - » K AUl 'i 
^frTE; P(l. ., ) , HT ( 1 ,)yS;.(3. l.CVClO'-) KAU03 


1 


F'0RM.AT( lX..3Eli3o / ) HLKOl 

13 FOP.^'IAT ( .1 HO .2 3H RESULTS i' '0.' i TKR 1 T ICM 9 1 F / ) 

14 Format (iHOal'rH Xl'^ GXLv F E vi/ ( 2X 9 1 w._ l :K. 6 ) ) 1 

15 Format (IHOkS^h ^...meral isi^c c i sps^-.ceflft vector >/( 2X9 lOt i i <,6 ) ) S 

16 FoR'-'AT UH 9 1.0X J.I fF>I- : r-U' CTIO. m' v^luE »l20- 0 /SX9l6'<I FUfjO'^IOTJ valuf : 

1 » 2 X,E 2090 ) 

17 For AT ( IHO 92 X 9 l<^ti ORTFk.^*^'ri.i/\LITY T L oT v E 19« B » 5X 9 . I 

: 1 6 H ’■■. V L c,\' G 1 M 9 1 9 .. !■.>. > B X 9 1 Vi >'’ 0 '/ c V EC'i 0 R ♦ E 1 9 „ n) . 

\H FOPsinT ( IrlO 9 ] 2 h :y\OVb VECTCA, s / ( . 2 X 9 1 0 L 1 1 « ^. ) ) I 


19 Fop;‘AT (1H094‘^I1 GnrjERALlSSO m I SPL ACEbt; -T VEC fOIM ST AClT I l-^G POJMT ) » 

1 /(2X.,1 OF-lloA) > 

20 For .aT (.lFl'>2X9;f2H SEAR^^ti FO.R E I GLi‘‘.Vr:CTOR TKK , 1 2 » ?X #3 iH ■ AN.C 

icoraiespOi'IOing eigenvalue 9 //? I 

21 Format (/ 922 H test for vfxtcr space) 

?.2- Format .(/pIri-i test .for ciTADiEfij) 

3 For"' AT (/ 92 IH TEST FOR. •■■rOV.E VECTOR)’ ' • 

.eKECK fCP EI.UEMVALy;E , On* <. I 2 9/ (.i!.X>.10Fl 1 * 6) ) ' 

;2H ■ P«:$T,ART60 ). 

,4H iHHC-f^^rsrCtP'tST ..DL,$CEl^X F,OV!:: ■ 

;.X,39-H! ‘F%L:$.F| FOR .. CrGEMVTXira;^ .hOo • 9i292X9.31- 

ipORO-f.NvG , ■eiGEWAtU'!E,.t» /7; 9U !' • . ■ , ■ ‘ ■ ' . 


29 FOB^ilAT (///^ -IHO, i6il, YiUjUi:iS (2X, : 

30 ; X'OxfiyUr (///,_ I 6 li Xf^"4Ui:^riCX'i;3- .vrt-, /, (XX^ r-;. 
131 FORf^AT(2X,10[:il„^) 

152 format (2X, 20 I 5) 

00 2 22 I = l.-.MnF5 
DO 2 22 J = 1‘,MVFC 
222 l;x(I5J) = OuO 

00 3 33 1 = 1 -''VFC 
DO 3 33 J = l,:jVFC 
333 XM3'.( 1 O) =0oD 

SDlVr=loOE+00 

lTM0=2-}f-:'.'DFS 
AO COMtTM'JF 
VFC=1 

WRITE (6,, 20) VFC 
DO >'4V I^lyMOFS 
AA9 X.{ I:j=vFco(VFC,I ) 

xx^r«c 

FU ’ i '‘-0 = 0 

Do 901 I = lv."DFO 

IF ( ( A33 { X ( 1 ) ) ) « Ll- » XX ) 00 TO 901 

XX=X(I) 

F;I^:o^I 

9vl CONTI^JUF 

DO AA7 I=1,^'■0FS 
4A7 X( I) = X(T)/XX 

WRITE (.6 ,lf^) (X ( I ) ^ 1 = 1 silRFf. ) 

453 . ]tP(VeC.,EQ« 1) 00 TO 49 

CALL T'OSTS(i\OFSi,^;\/LC»!'.V/.CCiX.;LiX»XOX ) 

WRITE ( 6 j 2 1) 

WRITC (XBX( 1 ) 9 1 = 1 V JKJ) 

49 ITFF=0 
3jC CONTINUE 
IT = 3. 

NI TFR^lO 

CALL F^MDOi iMDFS9SFC>X, oTIJUF »Fl9G) 

NUMFP 1. = NUMF.R. 

■|l = ori'.o, . 

IK (VLC 1) 70 8 5 

caI.l tests (, nFS4:-VEc,r:VFcc:-,09. X9X(3X) 

VJRTTE ( 6922 ) ■ 

hrite (6 93 ;) (XBx( n 9 I =1 9 JX.O) 

4449 IF (VL-^: oMf. .\|DFS) 00 TO 
VALU(VFC)=^F1 
' OMLOaI vrC)=SOR‘T (FI) 

. Do 67 T = l9'^’DF$ ■ 

67 v:6:.cmVT;c;9T)=x.( r ) , . 

-■ ■( 6-914) (•K{ n9r=^9'5WS).9f.O! i);9l-l9.^3D'F& 

GO TO' 

.liS 'DO ,T="l9T‘)D;FS -. 

5’0' 

.. WRLTE' ihrii)' ITER 

2c>a CACE ■FA•lXD;G^(KpF;S,9•^5•f^ . ' 


60 


NUivi£)R2=;FaMER 

0)\i2=3Dj:/i.- OM 
KU'^l-OoO 
DFM«0,'^ 

DO 51 T = l»'SnFS 
iMUM=MU'*+X( I ) ( I ) 

51 D^m = 0I:>i+X ( r ) ^(-63 { I ) 

. U = (iMUy.rF2^fDE'')-(f'‘U’3f-DtiVC' 2 .) 

V=: ( f:u, ‘Fr^-»^DC.MC-/!1 ) ~ ( i'-UJ-'FK I* Dr ■ ■n‘-;2 ) 

W= ( MU^-^fDEPO :i ) - (rU -iE-^l^DF'M 
V = V/U. 
b.'=l«7U 
u = un 

IPHTbRcGT.^l vAMDf ITCR^LToNITHR) GO TO 120C 
;■ I'/R ! T,E { 69 I) UyVj '»i 
iZ JO .TS1 = -V/ (ZoO-'t-U) 

. T<i2=( 5) / ( 2oO"'^'J) 

TS3«-TS2 " 

ALP,HAI = T31+TS2 
AUPMAi?'-TSl + fs3 

.l.F(ijtRvGT.<,,loANO«lTLR0.LT„{\!lTtPv) GO TO ?20C 
■. 'W;R,itfe,' ) A'LPHAl^ALPHAa 

22i;0' ip I ( AUPHAI vLTo 0^0) oA^^Uo (ALPHA2 <.LT» OdOl) GO TO 76 
IF ({ alpha loGTfiOdO) eAwDu (ALPilAicoTc 0-0) ) GO TO 60 
ALPI iA = AMAXl ( ALPHA1,ALPHA2 ) 

GO TO T'' 

60 ALPl!A=A'^lMl<ALr'MAl>ALPHA2) 

7( .IF(ITc^^GTdl.dANu«.ITt:RoLT«NlTi:p) GO TO 3300 
WRiT.E(6>l) ALPHA 
33.00 GO TO 71. 

76' DO 7-9 

79 XekI I. ) =X( U+ALPHAl^f-Sl I ) 

Call fapdg (noi-$,sp»xli-s^fal5Fo,fi*g) 

DO 60 T=1»H0FS 
bO XFi'\l n=XU )+ALPHA2*St n 

Call Fa 'DG <hiDF5 9SRs>XLlM9'..FAL3C-.»F2>C) 

IF( IT3r^Ct<1 uARl;.,ITER. LT.,f'ITER) GO TO 4200 
v^RITF (69I) F1,F2 
A2- 0 IF (FI r;LTo F2) GO TO 77 
ALPHA=;«LPHA2* 

GO TO 73 
77 ALPHAs^ALPrlAl 

IF ( ITLr*uTo 1 vKi'j-Ow ITlR«LTvIvITI^R ) GO TO '>200 
76 WrITF; {6-1-)- alpha 

52'jo do .-a i = i?'';dfs 

8.1 XFM( r).s^xn )+ALPhA-Jf>3( I )• 

0.0 3?,. T==l.>^«0(-5 

32 X.(T)=-XFM(n , ' . • 

( 6-» 25.) 

• S.T5 prof S. ?RV ii NVf-:CC 9 X XBX •) • ' 

6'#2.1).- • • . . ■ ;■•••' 

6.^1) ixBxt r) >TFa».3XJ )• ■ •• . • ' 


TA=*T;'\+G( I )^rG( I ) 

32 XtM( I ) =X ( I )+ALPHA^tS ( I ) 

Ri'-AX = OoC 
IMAX=0 

Do A-5]. I = l.>MOFS 

IF M^.BS(XFM( I ) ) ) cLtc. R-lAX) GO TO 431 

rmax=abs(xp:m I ) ) 

i:‘-iAX=I 

431 COVTIMUF 

IF(RXaX oLE. 10„0) go to 434 
NiNC='-iri 
rilM=IiiAX 
DO 452 I=ioNDFS 
452 X(,I )=XEi‘iTI) 

XX = X(F.tN) 

DO 450 I = 1:>MDFS 
450 XT I)-X( I )/XX 

Write (6,26,) 

,60 TO^ 

. 4.5.4 Caul, FAMt>G(NDFS,,5R,XEM,oTRUto ^FXETsGX^I-) 

‘ (vlUT'FRl^MUMER 
OFMCnlsDEFOr! 

1F( ITERoGTi^l^.AMD. ITERoLToFTTER ) GO TO 56 
IF(VEc'oEQul) GO TO 66 

call tests ( NDFS»NVEC5GVECCsXE.'1'>GX,XDX ) 

■ WRITE (6,2.1) 

write (691) ( XbX( 1 ) 9 I =1 9 JKJ) 

6666 CALL TESTS ( iMDFS » RVEC , NVECC , FXL|/! sOX v XBX ) 

WRitE (6 ,j22) 

ViRITE (Sal) (XBX ( I ) 9 1 = 1 9 JKJ) 

-86 GS = O,0 


AG-OoO ■ 

BE=^«0 

00 53 1 = 1 ,^'DFS 
GS = ns+S ( T ) ^GXEFi ( I ) 

AG = AG + F, ( I )-»!-S ( I ) 

53 6E = !-‘.E + GXLti( I )^tGXEKi( I ) 
BETA=:fc3E/TA 


AG=SQRT ( AG) 
UE = SQRT ( BE ). 


GO 


6100 


TEST = G5/( AG'>tBF) 

I F ( I T Ep? o GT 0 1 O' AND « I T ER « L T » N I TER ) 

WPsI'TE. (6 9 17.) TEST^LEtAG 
6aC0- DO 54 I^UNDFS 

<14 S( I ).- = -GXEM(- I l+iiEtA-X-SI I.) 

LE (VFCrt EO-I ) GO TO 620;0 , ^ 

I F ■( I TEr « GT (. 1 <> AND ^ ITER« Lto K I TlH ) GO TO 3^00 
GALL TESTS ( N'Dh S vriVEC a^VECC , 5 BX » XbX ) 

..(S?2:3 ) . ., • 

' -WRITE' ■{,6'9'i). (.XBK(.ri ?r = l»:JKJ)' 

. .v:^ 0.0 i F- ( I T liR GT .« 1 a I T E i?. «.:L:T p .M I T E U ) GO . .TO 

■WRlTETSaTS) 


B20Q 


DlF=Fl-FXe''' 

IF {Ilf uLT» GO TO ^rO- 

‘-•6 OP_l = :>IF-"<'5DTFF 
OFU = DtL/10..0 

IF (DEL vLTo FI) GO TO -'iCO 

IF ( (OFb LTu FI) ( /-Dcf { T.-ST ) ...T j a)) Gu TO 

99 DO 55 I = ljr-:OFS 
G( I ) =Gxr:u T ) 

;5 X( I )=XFM( I ) 

FI^FXEm 

ITEF;=^TtFR+1 

IF (ITFF oGT. IT^fD) GO TO 0 0^ 

IF{ I TFRcGT.rUTER ) IT = IT + 1 
L'rTER = 10^HT 

1F(1T'Ii.Od(jTi»J, i;-AMiv'L,ITc.RKjLl>j’\‘ITL.R) GO TO 
WRITE (6 ?13) ; ITER 

go: to 20g 
3GG Gr.MTT'UlE 

'Write (6 >14) (xeg(.i ) ? i = i,i-iDFG) ? (gxff'( 1 1 ? i^i ^ .dfs) ^fxe 
WRITE (6 >27) ,, VEC , 

G<5 TO' ^S'O 

ivO COf'lTtwuF. 

Write (6>14) (XETK I ) > I = 1/'‘CFG) s (GXEI’ ( I ) U = 1 sIvDFG ) gr-XE 
VALU(VF:C)=FXEf"; 

O.EGA (\/EC) =sa'U ( FXF'-I) 

DO 4.C1 I=1vI-!DFS 
^01 ytCO(VEC > I )=XEM( I ) 

Call }<x/tX.(MiDFs>SR ^xem, ..falsEo ) 

DO 6-tt ■ 

.$0 CV( l)^A/V(I )/r13( I ) 

WRIT-E(-,A^24) VEC > ( C\/( J ) » J = 1 si'/'ji-G } 

^.:.a iF- (VEC «F0<. r-VFO GO TO 69 
VFC=\/-r+i 
V<RTTF (6 >20) VEC 
^•■CU CC-.‘TjMJF 

CALL KX''X(i'IDFS>SR>XE!">oTRUE. ) 

JJ=VEC-1 

DO 5 8 I = 1>--)DFS 

uX( 1 >JJ)^F,l( 1 ) 

DO 40 T = 1>^>PF3 

40 F( 1 )=VEC 0 ( vr.c,i ) 

IF (V C .GT- 2) CC TO 57 
DEtn^O.,-^ 

Do 64 I = 1>'-'DFS 
6-4 DEM1-DE'‘'I 1 + P { I ) 'T-jX 1 1 > 1 ) 

DO' 65 

^ 65 X ( I.)=P( f ( I > U 

GO TO 92L ' 

I 57 I.F I VEC u-GT'.s . 3 ) . CO T-O 2' .<9 ■ ■ ■. ' 

■ D:Ec-a=- 

D0 ?9 ] rN'DFS ' 

^ FO- DE’;41=ii^'C 'i'i+'BXI I ».l ) (.1 


65 



Xf N(l9l )=lo ■ 

/■'AS 


Xf^L{i ,;m=-::L 

M/*-;A5 


Xf- N (29]) =-02 



XNN ( 2 9 2 ) = 1 . L 

/f^A5. ' 


GO TO 72 


22V 

VEC=VFC-1 



Call T'MCL‘;'(,’ 

L FS 9Ly XNr ,CF ) 


VFC = Vi-C + l 


72 

CALL FROD(i.;D 

FSs. vrc^p 96'', f X,. V X7.-.1 , , 


Do 61 I=1»0DF5 

61 

X( J )=P( T )-<* , 

^'P( I ) 

92 

xx=o.,^ 






b(3 9v2 1 = 1 v-: 

OFS 


IF ( ( AbS (X( I 

))) Ltr-ci XX) C-.^ TO P.,2 


XX'^XII), 



.•vliv=[, , 


9L' 2 

CGfJTINUF , 



PO ,62 T'=1>^'DFS 

■ ;62 

X( I)=X( I) /XX 


,i' 

WR ITE ;(6 ?1";) 

(X ( 1)91=1,, 'DF3 ) 


JKa=VFr,-i 



CALL TESTS ( 

'■;! D F S » V ;.: C > 0 v/ E C C u X 7 X I.1 X ) 


WRITF {6,211 



• \r-: I T ( 6 » 1 ) f x.'ix ( I ) v I - 1 jK j ) 

‘il-.'O 60 TO 

o:^ wZ-Mte. (6 92f') ( Vf-.LU { I ) J I 1 9'MV:.C ) 

■ v/t-x; I TE (61-30 ( OnLOA ( I ) , I = 1. , W CC ) 

end 


siiDFTr r/'i-nc 

subP.oui 1 r :: (-■.‘■■. ''bo ( mof:-,', sm , ? lt g -> i- , c: i 
INTFGFP SR.Vf-Cf 7bF 
KllAL ^'U-iER^'SiJr j -G 'MG^, ' ■'■ G 
L0GIC.''‘L LFG 

Co- Oi ■ /SI N G / Fi 9 P I.'. I- 9 L b I- 
Cc /PATHAIC/MVCC 
0 ‘'Vi-t LOGIC/ AA>rjV. 

COL* -G-'/CDLOCL/i.UMGP 
Co:- O-'-'/OBLnCK/Vt'C 
CCi- ON* /C. BLOCK / ax 

. CC'-- 0'^7btnXVCK/-:[M9 ,T -:C9lT£TR , 

GCf'- QP'/HARI/IFF 

D I ; fEnS lC-‘. X ( 1 O*' ) 9 G ( 1 b- J ) -i A A { 1 '* 0 ) 9 ha {• X ( 1 C b » 1 0 ) » m { too 9 1.0 ) -7 

] O y 1,2 ) V. B X ( 1 0 V 9 1 0 )7 h'> -H :>) G ( 1 0 ■; ) -i N 6 (. 1 0. • - -"C ( ) ■ ' • ' 

ut; ■( 6t:irv.,p; j ■' 



f; 


• iCAU*)0 
KAUO^ 


f 

•i 





65 


DO 33 T=l5NDFS 
33 N( I 9K)=BX{ J ,, J) 

31 CONTIMUF 

vec=ll+i 

GO TO 11 
12 Jj=VFC-l 

DO 26 T = l9Fr)FS 
26 BB ( I ) :r 6X ( i JJ 5 

gall In gen t NDFSsNWF.^Bb) 

DO 27 I = 1>NL';FS 
27 N { I ,VEC ) =BX { I j JJ ) 

11 Continue 

CALL Prod CNDFS 9 NVFC 5 ' J « 5 UN 9 NO . f 'f' NC . ' N NNO; J 
Do 37 I=l9MDrS 

37 G( n=G( I )-MNN\iG( I ) 

38 G (MIN) =0.0 
■ 2 A' RFT'UR'f’'' 

END ' ' 


$IBFTr PROD - ^ 

Subroutine p r o d ( n i.'j h s j . v l c ? * n > i \ n 5 r*. o ^ ■ i''! a 1 1 ^ g ) 

THI-S INTEFINAL SU-ROUT I f'X PERFOR'-iS TMi:. rMi^EL r-ATRlX VECTO,, MULTlP'-i^ 
CATIONS {r'(NTN}“l INT >G 
real: D 9 H N a N G 9 M N N G , N N ' ! G 
IMTEgEp? vec 
•c OMMO ^7 O'B L-O c K / V F C 

Dimension o( 100 ) jN c loi' v in > &( 10:; ) ? .^ang i ico ) oNN:\ng ( 100 ) jNn ( 100 ? 10 iicauo 

■Do 31 I=1,VEC 
SuM=^y n 

Do 32 J=19NDFS 

32 SUM = SUM+N { J 9 1 ) G ( J ) 

NG(I)=SUM 

31 Continue 

DO 33 T=i>VEC • , ■ • 

SUM=0d^ 

Do 3A J^UVEC 
3A 5UM = SU5A+NM (I ^ J ) ■'•''I'lG ( J ) 

M NHr-( I ) =sU'Vi 

33 CONTTNijE 
DO 35 I = l,aNDFS 
5UM=0o^ 

. DO 36 J=laVFG 
36 SUM~SUM + N( J 9 J)*NNNG( J) 

N'MNM6M)=5UM 
35 CONTINUE 

RETURN , 

■ BND' . . . ■ ' ■ . • 


r. o r, r, 


‘•>IBFTC KX''.X 

sub^iourr^E :;x:;x( .nrs^iir; 

5UU\..)UTI..'t: H'.bT UL 
p>Y Ti-./ r.[.t l.-i 

U)X Ai-J ( )X VtCTv' O* 
I'']TJ;G^:p Sr^ROp 

X^s.'T'v L 

Co*-”- 0^7s^5 1'’iCl-i/\7;'0p,L[;p' 
IS r.f p(_;- v 

Is XU'.JCt o(- 

1 5 i^U’bLO 01- L.^ ,oii 

L0G1CA|_ Lbr'- 
C'OM' Or /BoLOGrv-vAj L-b 
COM-’Oh’/HAR I / iff 
L),Il EM sr^ '.i (,:.2 9 12 ) 

12 ) 9 r.s2 ( 10;;'9l2 ) 9 AA( J • •; 
ASSIGM 5VSTh:.v. t^^wfo-LS 01 
lFF=r,lFFH:l 
, ' R FF/ TND’ ,3 '; 

^ ';ia' Tpl.v’x 

C>Q. ;i;,l:, 7=i,,;|:rOFv, ' ; 

Y ( i 9 J ) r. ■'( i» '"J 

JJ = SR ( T 9 J ) 

IF (jj} 1109 n»n.'- 

‘ILO Y(IU)^X(JJ) 

11 CO-iTpHUF: 

10 COMTIMUE' 

: Calculate <k)x «M-i* ( 

Do 12 .rsl,,M 
I F( IF F'^GT* 1 ) GO TO .f;o 
Call ) 

Write (3) ( (kpij^x) oi=i9i 

GO TO 27 

>A READO) ( ?'sp^lsi 

>7 Do 13 J=1»RDF 

suHi= ';*,o 

Do lA K.= l»nCF 
su^u= SUv.]+K'Mw,X )rr( I 9 
L4 SUM 2 = SUN 2+y, ■; ( J 9 X,) Y t I 

TSl(l9j)=. SUMl 
3 TSPdfJ.)*’ StV,2 

12 CONTI.MUF 

■Do 3 a ip.l^rdFS 
aa( I )”-:7 0 

13 Bi3(i);=3,0 

DO' 19 Kp1*x. 

DO 20 K.vfei'i,.RvJF' . 

I — SR ( cs ii“s K ) 

IF IJ) 2.1., 20, .M 

AA( I l) + Ts.(. t'‘ ,<!<:,) ' ' 

O-d '4'C‘SL 

J.r'lJ's , 

CO^^;TTf.'vU:E 
IF. (:«'i(jnT> :V 


'3 R A S 1 = r? AS 1 + R ( I ) -X- c B ( 1) 

^0 24 T=l 9 MDFS 
DB ( I ) =BD { I ) / ( ( |<AS1 ) *->5-0 3 ) 
^2 KeTUR^• 
t:ND 


■3 IE FTC INGEM 

Subroutine ingemc-^dfs^a^;-;, >x; 

^ ^UI5 iN'TERf-iAL SUbRC^UT I ,\'E GENER^hTES ThE TmvppS>^ 

C NVEC Is number of 1 1 GiiX'VALUES-t I TENVrCTORF Drr.iirD(j: 

integer VEC 

real f' ' N V 9 M G 9 K> N F G 3 M G 

Co>.h.'ON/PATHAK/NVEC 

comhon/dbuock/vec 

DlijvlENSlOEi PK{ 10,0 ) , T i ( 100 ? i ? J oAl ( i or. » i a ) , ^2 ( ) » X ( 

■ ■■ . N' (10 0 >' 1 O' ) 9 N 6 ( 1 O' e } ■ 9 H H G ( I O O ) > N N N rw l n n j 

■ CALL' P.ROLi ('NDES y.NVFC ? X s N.» Rjv » .’■■G s FMfVG 9 i.M.MrVr; ) 

bo 10 i = i 5 N’Drs 

10 PM( I 5 =X{ I )-FfNWG( n 

Ac=Ov'‘ 

DO 11 J = 1 >FB-FS 

11 AO^AO*i-pM ( t ) ( 1 ) 

AO^l.O/AO 

K-KK-^^VilC-l 

DO IZ J^-l^KKK- 
90 13 

'13 rM{ I , T )->^'v''^'f'‘G( J J 

3 <1 »^i^lMTTMnf 7 . 

On i/«. 

na IF I - 1 j 1 /'''.'^ 

15 Ai / T , I ^ ( I f, J) + AC^:-T":( I , J ) 

Camtt».!i.“ 

: l^ltKKK 

1 ^ ^ J =_ * C 4 ^ 

no 17 .1"], »KFK 
if .v.\!( I » J)=/.l‘( I V 
JJ=VEC 

DO 18 T=l9RKK 
MN( JJ> T ) =A? ( T i 
18 NWl; T-,,l L)^AO( > J 
NW( = V-? 


. 
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n rr T) I*"* M 
crki r» 


r 

r-!; 1 R F T r T F c T c 

Si IRPnt IT T MF TCCT C ( Mrxc c ^ j V > 

ImtffFo \/Fr 

, CnMMnM/nR( nrt/ /\/Fr 

‘ DTMFMCTOM Y ( t Rn \ jPY ( 1 nr* ^ 1 , vnvy { t j 

C0W1=;A^A 

nr»l 0 T.= 1 j, M'r^ir c 

10 tOMl =rnM1 +v ( T f T ) 

CF\M1 =1 ( t FAKII ) 

■ nni 1 j-1 J J J 
' , Si = A ^ 

■ ■ ' ■ on 1 7T z; 1 J 

19 Sum- fmmxfv ( T J t ) J!.v ( I J 
Xd V { J \ -Cl I^jt 
11 rAMTTMMT 

HAI Q I-*- jJJ 

y°y ( J )=:yp.y. ( J )/cc?\a 

Q r Tl I Ci » t 





